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Abstract 

We extend the technique of constructive expansions to compute the 
connected functions of matrix models in a uniform way as the size of 
the matrix increases. This provides the main missing ingredient for a 
non-perturbative construction of the (f>% 4 field theory on the Moyal four 
dimensional space. 

1 Introduction 

Constructive field theory build functions whose Taylor expansion is perturba- 
tive field theory |2] . Any formal power series being asymptotic to infinitely 
many smooth functions, perturbative field theory alone does not provide any 
well defined mathematical recipe to compute to arbitrary accuracy any physical 
quantity, so in a deep sense it is no theory at all. 

In field theory infinite volume quantities are expressed by connected func- 
tions. One main advantage of perturbative field theory is that connected func- 
tions are simply the sum of the connected Feynman graphs. But the expansion 
diverges because there are too many such graphs. 

In fact connectedness does not require the full knowledge of a Feynman graph 
(with all its loop structure) but only the (classical) notion of a spanning tree 
in it. To summarize constructive theory, let's say that it is all about working 
as much as possible with the trees only, and rcsumming or hiding most of the 
quantum loops. This is the constructive golden rule: 

"Thou shall not know all the loops, or thou shall diverge!" 

However the constructive program launched by A. Wightman and pursued 
by J. Glimm, A. Jaffe and followers in the 70's was a partial failure because 
no natural four dimensional field theory could be identified and fully built. 
This is because only non-Abelian gauge theories are asymptotically free in the 
ultraviolet limit. But ultraviolet asymptotic freedom also means infrared slavery, 
and non-perturbative long range effects such as quark confinement are not fully 
understood until now, even at a non-rigorous level. The constructive program 



went on, but mostly as a set of rigorous techniques applied to many different 
areas of mathematical physics [3J H] . 

Recently quantum field theory on non-commutative space has been shown 
rcnormalizable. The simplest such theory is the theory on the Moyal space, 
hereafter called (p^ 4 . Grosse and Wulkenhaar 5 overcame the main obstacle 
to renormalizability, namely the ultraviolet/infrared mixing, through the use 
of a new symmetry called Langmann-Szabo duality [6 . Following their initial 
breakthrough, a flurry of papers has appeared to extend this result to other 
models and to generalize to the Moyal context many useful tools and techniques 
of ordinary perturbative field theory. For recent reviews, see [7J[S]. 

It now appears that four dimensional non commutative field theories are bet- 
ter behaved than their commutative counterparts. In particular ^J 4 , in contrast 
to its commutative counterpart, is asymptotically safe [9l 1101 ITTj : the flow be- 
tween the bare and the renormalized coupling constant is bounded. In fact the 
graphs responsible for the flow of the coupling constant compensate exactly at 
any order with those responsible for the wave function renormalization. This is 
an exciting discovery: LS symmetry may play a role similar to supersymmetry 
in taming ultraviolet flows. 

Asymptotic safeness is in a sense much simpler than asymptotic freedom, 
and </>4 now stands out as an obvious candidate for a four dimensional con- 
structive field theory without unnatural cutoffs (although on the unexpected 
Moyal space) . 

But after |llj one main difficulty remained unsolved on the road to con- 
structive </>4 4 . Current cluster expansions used in standard bosonic constructive 
theory [2] are unsuited to treat matrix models with large number of components. 
To explain why, let us compare the large N vector <f> A model and the large N 
matrix[12 </> 4 model. In both cases the coupling scales as 1/N for a non trivial 
limit as N gets large: at order n in a graph there is indeed in both cases at most 
about n loops of indices. But in the first case the field has N vector compo- 
nents, and at a given vertex only two different colors can meet. Knowing only a 
spanning tree in the graph, it is still possible to sum over all indices at the right 
cost. To do this, at any leaf of the tree one can sum over the index which does 
not go towards the root and keep the other one for the next step. Iterating from 
leaves of the tree towards the root there is only one index summed per vertex, 
(except at the root, where in the case of a vacuum graph there are two indices 
to sum, leading to the final global N factor of vacuum graphs) . This procedure 
does not violate the constructive golden rule, as no loops need to be known. 

But a matrix model is very different. The field has N 2 components and at 
a given vertex four different indices meet. The scaling of the vertex is still only 
1 /N, but this is because each propagator identifies two matrix field indices with 
two others, rather than one. Therefore matrix models apparently clash with 
the constructive golden rule. The knowledge of the full loop structure of the 
graph, not only of a tree, seemed until now necessary to recover the correct 
power counting, for instance a single global N 2 factor for vacuum graphs. 

Since (f>^ is a quasi-matrix model with a large number of components in the 
ultraviolet limit [13 it is plagued with this constructive matrix difficulty, hence 
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seems unsuited at first sight for a constructive analysis. The difficulty persists 
in the direct space version [2] of the model, but in a different guise. In that 
representation, it is the non-locality of the vertex in x space which is impossible 
to treat with standard constructive methods, such as ordinary cluster and Mayer 
expansions with respect to lattices of cubes. 

In short a new kind of expansion based on a new idea is required for con- 
structive 4>\ A . This is what we provide in this paper. 

The idea is in fact quite simple. Matrix models can be decomposed with 
respect to an intermediate matrix field. Integrating over the initial field leads 
in a standard way to a perfect gas of so called loop vertices for this intermediate 
field. One can then perform the tree expansion directly on these loop vertices. 
All indices loops then appear as the correct number of traces of products of 
interpolated resolvents, which can be bounded because of the anti-Hermitian 
character of the intermediate field insertions. 

We take as an example the construction of the connected functions of a 
matrix model perturbed by a ~Tr^>*</><^*0 interaction. We prove as a typical 
result Borel summability in A of the normalization and of the connected 2p point 
functions uniformly in the size of the matrix Q 

In a companion paper |151 we explore the consequences of this idea in the 
more traditional context of commutative constructive field theory. 

Recall that it is possible to rearrange Fcrmionic perturbation theory in a 
convergent expansion order by order by grouping together pieces of Fcynman 
graphs which share a common tree |16L I17j . But bosonic constructive theory 
cannot be simply rearranged in such a convergent way order by order, because 
all graphs at a given order have the same sign. Resummation of the perturba- 
tion theory (which occurs only e.g. in the Borel sense) must take place between 
infinite families of graphs (or subparts of graphs) of different orders. To ex- 
plicitly identify these families seemed until now almost impossible. Cluster and 
Mayer expansions perform this task but in a very complicated and indirect way, 
through an intermediate discretization of space into a lattice of cubes which 
seems ad hoc for what is after all a rotation invariant problem. 

In fact the cluster expansion between loop vertices, although found in the 
context of matrix models, can identify such families also in the ordinary commu- 
tative case [15 . This simplifies traditional bosonic constructive theory, avoiding 
any need for cluster and Mayer expansions. We should bring in this way Bosonic 
constructions almost to the same level of simplicity than the Fermionic ones and 
explore the consequences in future publications. 

2 Matrix Model with Quartic Interaction 

The simplest </> 4 matrix model is a Gaussian independent identically distributed 
measure on N by N real or complex matrices perturbed by a positive -^Tr(p*(j)(f>*(f> 
interaction. The N — > oo limit is given by planar graphs. It can be studied 

1 Non- uniform Borel summability, taking A smaller and smaller as N — > oo is trivial and 
would completely miss the difficulty. 
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through various methods such as orthogonal polynomials [TSJQIJ], supersymmet- 
ric saddle point analysis [201 EZD H2 and so on. However none of these methods 
seems exactly suited to constructive results such as Theorem |3 . 1 1 below . 

Consider the complex case (the real case being similar). The normalized 
interacting measure is 

where 

dfi = Y[dm ij d%$ ij (2) 

is the normalized Gaussian measure with covariance 

< >=< $ij$kl >= 0, < $ij$kl >= SikSji. (3) 

For the moment assume the coupling A to be real positive and small. We de- 
compose the $ functional integral according to an intermediate Hcrmitian field 
a acting either on the right or on the left index. For instance the normalization 
Z(X,N) can be written as: 

Z(X, N) = J d^GUE^e-^ 10 ^ 1 ® 1 ^^ 1 ®^ (4) 

where dfiQUE is the standard Gaussian measure on an Hermitian field a R , that 
is the measure with covariance < cr^a^ >= SuSjk- The e _Trlog represents the 
Gaussian integration over hence a big ./V 2 by N 2 determinant. It is convenient 
to view M. N as M^^R-^. For instance the operator H = \J~j^[t®v R ] transforms 

the vector e m ® e„ into ^J^e m <S> X)fc cr fcL e fe- Remark that this is an Hermitian 
operator because a R is Hermitian. 

By duality of the matrix vertex, there is an exactly similar formula but with 
a left Hermitian field a L acting on the left index, and with [a L (£> 1] replacing 
[1 <S> a ]. From now on we work only with the right field and drop the R 
superscript for simplicity. 

We want to compute e.g. the normalization Z(X,N), which is the (Borel) 
sum of all connected vacuum graphs. We define the loop vertex V by 

V=-Trlog(l®l + l®iff), (5) 

and expand the exponential as -^y. To compute the connected graphs we 
give a (fictitious) index v — 1, ... , n to all the a fields of a given loop vertex 
V v . At any order n the functional integral over dv(a) is obviously also equal 
to the same integral but with a Gaussian measure dv({a v }) with degenerate 
covariance < cr^-cr^ >= SuSjk- We apply then the forest formula of [23] to test 
connexity between the loop vertices from 1 to n. The logarithm of the partition 
function or pressure is then given by the corresponding tree formula exactly like 
in the Fermionic case [17 . 
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Theorem 2.1. 

io g z(A,iv) = f;^{n[/ 1 ^ E ]}/<MK}»-M) 

n=l T ^fGT ^° ie,je,k e ,k 

where 

• eac/i fee £ of the tree joins two different loop vertices V v ^ and V v ^ , 

• the sum is over trees over n vertices, which have therefore n — 1 lines, 

• the normalized Gaussian measure dz^r({cr„}, {w}) over the vector field cr v 
has covariance 

< a ij a ki >= 6u6 jk w T (v,v' ,{w}) 

where w T (v, v' , {w}) is 1 if v = v' , and the infimum of the we for I running 
over the unique path from v to v' in T if v ^ v' . This measure is well- 
defined because the matrix w T is positive. 

This is indeed the outcome of the tree formula of in this case. This 
formula is convergent for A small enough! 

Theorem 2.2. The series is absolutely convergent for A small enough. 

Proof Consider a vertex V v of coordination k v in the tree. Because the a field 
acts only on right indices, and left indices are conserved, there is a single global 
N factor for V v coming from the trace over the left index. We can then from 
now on essentially forget about the left indices except that they give a particular 
cyclic order on V v . See Figure[l]for a tree on four loop vertices, hence with three 
lines. 



Figure 1: A tree on four loop vertices 

We compute now the outcome of the k v derivatives Yli=i sTF acting on V = 
— Trlog(l + iH) which created this vertex. Fix an arbitrary root line £q in the 
tree T. There is a unique position i = 1 on the loop vertex from which a path 
in T goes to £ a , and the loop vertex factor V v after action of the derivatives is 

2—1 i—1 
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where the cyclic convention is k v + 1 = 1, and the operator C(i, i + 1; a v ) — 
(1 + iH(<T v ))~ 1 (ji,j i+ i) acts only on the right index (it is no longer a tensor 
product, since the left trace has been taken into account in the global N factor 
in front of V v ). 

To bound the integrals over all sums j e k( ifL (which by the way are only 
about right indices) we need now only a very simple lemma: 

Lemma 2.1. For any {w} and {cr v } we have the uniform bound 

HI E MiA . l W ,i W ])l[V v \<N* (8) 

Proof Since iH is anti-hermitian we have indeed ||(1 + < 1. The 

product over all vertices of the resolvents C(i, i + 1; a v ) together with all the 
sums J2i e j t k t i e exact ^y forms a big trace of 2(n — 1) operators which turns 
around the tree (see Figure [2|. This is the key point. This trace of an operator 
of norm smaller than 1 is bounded by N. 



Vertex 



Resolvent 



Figure 2: Turning around a tree with four vertices and three lines 

It remains now to collect the other factors. There is an N factor for each 
vertex of the tree and a \—iy/X/N\ factor for each half line of the tree. Collecting 
all the N factors we get therefore a a iV 2 global, n independent factor as should 
be the case for vacuum graphs in this matrix $ theory, times A™" 1 . □ 

We can now integrate the previous bound over the complicated measure dvr 
and over the {wi} parameters. But since our bound is independent of o~ v and 
w's, since the measure dv(a) is normalized, and each we integral runs from to 
1, the bound goes through. 

Finally by Cayley's theorem the sum over trees costs nl Y[ v ( k The 
n! cancels with the 1/n! and we remain with a geometric series bounded by 
N 2 Yln>i ^" 1 nence convergent for A < 1. 



This completes the proof of Theorem 2.2 □ 



3 Uniform Borel summability 

Rotating to complex A and Taylor expanding out a fixed number of vertices 
proves Borel summability in A uniformly in N. 
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Definition A family fj\f of functions is called Borel summable in A uniformly 
in N if 

• Each /at is analytic in an N independent disk Dr = {A|RcA _1 > 



Each /a? admits an asymptotic power series J^k a N.kX k (its Taylor series 
at the origin) hence: 



r-l 



MA) = 5>Ar, fe A fe + i?Ar, r (A) (9) 



fe=0 

such that the bound 



\Rn,vW\ < A N a r r\\X\ r (10) 

holds uniformly in r and A £ Dr, for some constant a > independent 
of N and constants An > which may depend on N. 



Then every fx is Borel summable [24], i.e. the power series J2k a N,kjj 
converges for \t\ < A. It defines a function Bj^(t) which has an analytic con- 
tinuation in the N independent strip S a — {t\ dist (t, M. + ) < -}. Each such 
function satisfies the bound 

\B N (t)\ < B N eH forteM+ (11) 

for some constants Bn > which may depend on N. Finally each f^ is repre- 
sented by the following absolutely convergent integral: 

1 f°° 

fN(X) = jj e-iB N (t)dt forAeCfl. (12) 

Theorem 3.1. The series for Z{\ 1 N) is uniformly Borel summable with respect 
to the slice index N . 

Proof It is easy to obtain uniform analyticity for 3t\ > and |A| small enough, 
a region which obviously contains a disk Dr. Indeed all one has to do is to 
reproduce the previous argument but adding that for H Hermitian, the operator 
(1 + ie^H)- 1 is bounded by 2 for |0| < n/A. Indeed if tt/4 < Argz < 3?r/4, we 
have |(1 + iz)~ 1 \ < V2. 

Then the uniform bounds (10) follow from expanding the product of resol- 
vents in Q up to order r — 2(rt — 1) in A. by an explicit Taylor formula with 
integral remainder followed by explicit Wick contractions. The sum over the con- 
tractions leads to the a r r\ factor in (10); in our case the constants An — K.N 2 
actually depend on N but this is allowed by our definition of uniform Borel 
summability. □ 
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4 Correlation Functions 



To obtain the connected functions with external legs we need to add resolvents 
to the initial loop vertices. A resolvent is an operator C(a r , mi, ^2), which 
can depend on only two indices because in a matrix model every entering index 
must go out. The connected functions S ,c (m 1 ,...,m2 P ) therefore depend only 
on 2m, not Am indices. They are obtained from the normalized functions by 
the standard procedure. We have the analog of formula [6] for these connected 
functions: 

Theorem 4.1. 

jxi^ ( _ r l ~| ! 

du T ({a v },{w}) 



^...^H^fWnt/ 1 ^ e ]}/- 

tv n=i t leT J ° u,jt,kt,h 

n [Si t i t s jtke - s v(i) v >u) ]\\n v v n Q^^r,^^))} 

tZT 5a iU Sa k e ,h J 1 " r=l J 



(13) 



wht 



• the sum over n runs over the pairings of the 2p external variables into 
pairs (^( r ,i),^( r ,2)), r= 1, 

• each line i of the tree joins two different loop vertices or resolvents V v ^ 
and V v i[t), 

• the sum is over trees joining the n + p loop vertices and resolvents, which 
have therefore n + p — 1 lines, 

• the measure dvr^a 1 "}, {o~ r }, {w}) over the vector fields {cr Q } has covari- 
ance 

< a ij a kl >= SiiS jk w T (a,a' ,{w}) 

where again for a, a' E {v},{r}, w T (a,a',{w}) is 1 if a — a', and the 
infimum of the W£ for £ running over the unique path from a to a' in T if 
a ^ a 1 . 

This expansion is convergent exactly as the initial one and we get: 



Theorem 4.2. The series (13) is absolutely convergent for A small enough, and 
we have: 

\S c (m 1 ,...,m 2p )\<K(2p)\\N 2 -P. (14) 



5 Further topics 

5.1 Symmetric or Hermitian matrix models 

Interacting GOE and GUE models can be treated along the same lines. Let 
us consider for instance the same model than (ITT) but with $ = $* now an 



8 



Hermitian matrix. We have no longer a canonical distinction between left and 
right indices so that the intermediate field operator acts on both sides, but it is 
still anti-Hermitian. The vertex operator (J5J is therefore replaced by 

K = -Trlog(l®l+|yA[ ff0 i + i 0(T ]) ; (15) 

so that each loop vertex is no longer simply proportional to N because of e.g. 
the left trace. But any tree is planar so one can still draw the tree between 
loop vertices on a plane, as in Figure 3. The total number of traces of products 
of (1 + iH)^ 1 operators for a tree on n vertices still remains n + 1 by Euler 
formula. Indeed Euler formula says 2 — 2g = V — L + F, where g is the genus 
and F is the number of faces, each costing N. But graphs of genus as those 
of Figure 3 contain 2(n — 1) vertices (of the cubic type), and two kinds of lines, 
the n—l lines of the tree and the ^2 k v = 2{n — 1) resolvent lines. Therefore 
F = 2 - 2(n — 1) + (n — 1 + 2(n - 1)) = n + 1 so that all the results of the 
previous sections remain valid. 




Figure 3: A tree with five loop vertices joined by four tree lines, eight cubic 
vertices, and eight resolvent lines which is a contribution in the Hermitian case. 



5.2 Genus expansion 

The genus expansion which lies at the root of matrix models can be generalized 
constructively. We can indeed expand the resolvents on the external loop which 
turns around the tree in Figures 1 or 3, and Wick-contract one at a time the a 
fields produced. If we were to perform this to all orders the expansion would 
diverge. However we can also contract until a fixed number of non-planar cross- 
ings are generated, and then stop. We call this expansion a "rosette expansion". 
It does not diverge and allows to extract the 1 /N expansion up to a fixed genus 
g, plus a smaller remainder. For instance for the normalization one obtains a 
representation 



Z(X,N) = N 2 



J2N- 2k a kA + R g+1 (X,N) 



k=0 



(16) 



where afc(A), the sum over all vacuum graphs of genus k, is an analytic func- 
tion of A and i? ff +i(A, N) is a convergent series whose sum is bounded by 
0(jV-2(s+i)) and 

is again Borel summablc in A uniformly in N. 
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This remark is essential to construct 4>\ through this method. We need 
indeed to identify the planar contributions with a single broken face and two or 
four external legs because they are the only ones which need to be renormalizcd, 
and also the only ones which can be renormalized (because only planar graphs 
with a single broken face look like Moyal products when seen from lower renor- 
malization group scales [5]). It is therefore essential to have a method which can 
extract them from the rest of the expansion without violating the constructive 
golden rule. This can be done through the rosette expansion sketched above. 

5.3 Decay of correlations in quasi-matrix models 

To fully construct §* 4 we have to take into account the fact that the propagator 
of (fy 4 in the matrix base does not exactly conserve matrix indices |13j , except 
at = 1, where is the Grosse-Wulkenhaar parameter. 

It is therefore essential to show not only uniform convergence but also decay 
of connected functions with respect to external matrix indices in this kind of 
models. This should not be too difficult using iterated resolvents bounds, as is 
shown in |15j in the case of ordinary cf> 4 on commutative space. 

5.4 Multiscale Analysis 

To fully construct (fi\ 4 we have also to generalize the single N analysis of this 
paper to a multiscale analysis such as the one of |25j . This requires to optimize 
as usually the tree expansion over all the scales so that connected functions of 
higher scales are always correctly connected through the tree. 

In fact the (p^ 4 can presumably also be built as easily in x space represen- 
tation by a slight modification of the matrix argument. Indeed a Moyal (j> 4 
vertex can be decomposed in terms of an intermediate ultralocal real field with 
a Tr<^> * <j)-ka interaction. This can again be done in two ways by duality. The 
new vertex is anti-hermitian again as a kernel between the </> and <f> points. The 
bosonic covariance of the </> field is a Mehlcr kernel that can be easily broken in 
square roots. We obtain loops of Mehler kernels sandwiched between operators 
of the (1 + iH) -1 type. We expect therefore all constructive aspects to be also 
doable in x-space [25] , 

Since our loop vertex expansion seems very well suited to treat both large 
N vector and large N matrix limits, we expect that it is the right tool to glue 
different regimes of the renormalization group governed respectively e.g. in 
the ultraviolet regime by a small coupling expansion and in the infrared by a 
"non-perturbative" large N expansion of vector or matrix type. This gluing 
problem occurs for the vector case in many different physical contexts, from 
mass generation of the two-dimensional Gross-Neveu |27j or non-linear a-model 
[28] to the BCS theory of supraconductivity [29] . Confinement itself could be 
a matrix version of the same gluing problem [12 . All such gluing problems 
have been considered until now too complicated in practice for a rigorous (i.e. 
constructive) analysis. We hope that this might change over the coming years. 
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